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Abstract 

Available experimental data on decay rate and polarization are used 
to investigate non-factorization contribution to processes of the kind 
B — > ATV') s-nd B K*'ip using five theoretical models for the formfac- 
tors. Using the knowledge on non-factorization gained from B decays 
we study the processes Bg — > {r]' ,r])ip, and Bg (j)tp where experimental 
data are very limited. 

(PACS numbers: 13.25.Hw, 14.40.Nd) 



1 Introduction 

The 1/Nc expansion of hadronic matrix elements has been a very important 
approach in the study of weak non-leptonic decays. The leading terms of this 
expansion are factorizable into simpler ones, whereas the next to leading terms 
are not completely factorizable. 

In the standard approach Fierz transformation and color algebra are used 
to transform the non- leading contribution into a factorizable part, which is 
added to the leading terms, and a non-factorizable part which is neglected [0. 
This is called the factorization approximation which is extensively used and 
sometimes works well. 

It was shown by Gourdin, Kamal and Pham that factorization approx- 
imation in all commonly used models of formfactors could not account for the 
longitudinal polarization in i? — >■ K*il), Pl/P, and the ratio {B Kip)/[B — *■ 
K*ip). Subsequently, it was shown in 0] that inclusion of non-factorized terms 
enabled one to understand both data in all the commonly used models of 
formfactors. 

Our aim in this work is to investigate the non-factorization contribution 
to processes of the form B — > Kip and B — > K*ip using the experimental data 



available on decay rates and polarization. This is a more thorough analysis 
than that presented in 0. Then, in the light of this investigation we study 
the processes — >■ rjip, —>■ ri'ip and B^ — ^ (pijj. 

The layout of this paper is as follows: Section 2 deals with the definitions 
and introduces the five theoretical models of formfactors we employ. Section 
3 contains an analysis of non-factorization in i? ^ Ktp and B K*iIj decays. 
Bg decays are discussed in Section 4. We end with a discussion in section 5. 



2 Aim 

In this work we are interested in nonleptonic B decays of the kind {B, Bg) — >■ 
PiV)^, where B is {B+ or , P is {K+,K°,r] or V) , V is Ik*+,K*^ 
or (f)) and ip represents ip or iIj{2S). These processes have similar flavor flow 
diagrams - see Fig. |l] - but different spectator quarks. 

2.1 Decay Rates and Polarization 

Using the effective Hamiltonian that contains the short distance Wilson coef- 
ficients Ci and C2, the decay amplitude for such processes is written as [§] 



A{B^P{V)^) = {P{V)mf\B) 
~ V2' 
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Nc is the number of colors and A" are the Gell-Mann matrices. 

At this point, the number of colors, Nc = 3, will be taken seriously. Also, 
instead of neglecting the non-factorizable terms ( second term in (0) and any 
non-factorized contribution to the first term ), we paramerize them as in [Q. 
The parametrization is done in such a way that we can conveniently combine 
the factorizable and the non-factorizable terms. Explicitly, this means writing 
the following: 
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where 



= iVB - VP(V))t, = {P^)t,- (9) 

The polarization vectors and t]^ correspond to the two vector mesons ip and 
V , respectively. In (|^), Cp arising from the mixing of r] and r]' (see reference 
is givin by 



{cosOp + ^smOpj P = 7] 
(^-^ cos 6p — sin 6 p^ P = rj' 
1 



(10) 



P = K+ 



where the mixing angle is taken to be 

Op = -20°. 



The non-factorized contribution to the first term in (|l|) is parametrized 
analogously to (0) and (§) by writing F^^^^^ , A^i^^^ , etc. in place of Fi^^^^ , 
Af)^^ etc. . 

Substituting (§ - ||) into the decay amplitude (Q), we can calculate in 
a straight forward manner decay rates for the processes B P{V)ip and 
polarization for the processes B Vip. The results are presented below: 
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Subscripts L and T in (0) and (|T5D stand for 'longitudinal' and 'transverse' 

and 00, ++ and represent the vector meson helicities. 

In (|ll| - |T5|) we have introduced the following dimensionless parameters: 
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Furthermore, x and y represent the following ratios, 
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2.2 Formfactors 



Before proceeding to calculate the decay rates and polarization we need the 
values of the formfactors. In this work we consider five theoretical models. 
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The first is the original Bauer-Stech-Wirbel model |^ (called BSW I 
here) where the formfactors are calculated at zero momentum transfer and 
extrapolated to the desired momentum transfer using a monopole form for all 
the formfactors Fi, Ai, A2, and V. 

The second is a modification of the above (called BSW II here). In this 
model the extrapolation from the zero momentum transfer is done using a 
monopole form for Ai and a dipole form for Fi, A2, and V f^. 

The third is the model of Casalbuoni et al. and Deandrea et al. ||^, 
where the normalization at = is obtained in a model that combines heavy 
quark symmetry with chiral symmetry for light degrees of freedom. We call 
this (CDDFGN) model. Here, all formfactors are extrapolated with monopole 
forms. 

The fourth is the work of Altomari and Wolfenstein (AW) in which 
the formfactors are evaluated at maximum momentum transfer based on a 
non-relativistic quark model. The formfactors, here, are extrapolated with 
monopole forms. 

The fifth is the non-relativistic quark model by Isgur, Scora, Grinstein and 
Wise (ISGW) I^Of. In this model the formfactors are calculated using harmonic 
oscillator wave functions for the particle states. This results in an exponential 
dependence in (momentum transfer)^ for all formfactors. 

The predicted formfactors in these five models relevant to the processes of 
interest are shown in Tables |^ and 



3 Analysis 

In this section we study the non-factorization contribution to the processes 
B — > Kip and B K*ip. This has been done in but our analysis is more 
thorough. First, we need to assign values to the CKM matrix elements, the 
Wilson coefficients and the decay constants which we take to be, [[TT], |^ 





= 0.974, 




= 0.04, 




= 1.12 ±0.01, 


C2 


= -0.27 ±0.03, 




= 0.384 ± 0.014 GeV, 


ftp(2S) 


= 0.282 ± 0.014 GeV. 



3.1 B ^ Decays 

The processes considered in this paper are color-suppressed. This implies 
that the first term in (|l|) is proportional to 02 which is rather small, 0.10 ± 
0.03. In contrast, the coefficient of the second term in (|lD being Ci is an 
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order of magnitude larger. Hence, as can be seen from (|TT]), even a 10% non- 
factorization in xfi causes the decay rate to increase four times, and with a 
20% non-factorization, the decay rate becomes nine times larger than the value 
predicted in the factorization approximation. 

Let us first consider the process i?^ K^ip whose decay rate is more 
precisely measured than of the neutral mode. The decay rate formula for this 
process, Eq. ([TTD , can be rearranged and written as 



- ^ 



{T{B+ K+ijj)} 



1/2 



(78.422 X 1012 GeV-' sec-i)V2 K^Kfeas^FiK ) 



(29) 



By substituting the experimental value for the decay rate ||lT| we can calcu 



late the non-factorization contribution as a function of the formfactor Fi(m^). 
This is shown in Fig. ^ where the curves represent xfi with its variance, 
Xfi ± o'xFi ■ The allowed region hes between the two curves. 

This graph shows that if a model predicts a small value for the formfac- 
tor, a high non-factorization contribution is needed in order to explain the 
experimental data. The five models considered in this work predict formfactor 
values in the range (0.55 - 0.84) ( see Table and the dots in Fig. ^ ) which - 
according to Fig. ^ - need about ( 10 - 25 % ) non-factorization. 

In a different manner to exhibit the same physics, we show in Fig. |^ the 
relationship between the branching ratio and the non-factorization contribu- 
tion xfi for the processes — > K~^ip, K'^ip, and B~^ K~^ip{2S) for 
the five theoretical models we have considered. The horizontal lines represent 
one standard deviation bounds on the branching ratios [0. 

Fig. H strongly suggests that factorization {xfi = 0) does not work, and 
that non-factorization contribution could help to explain experimental mea- 
surements. The first two graphs in Fig. |^, representing processes B~^ — > K~^il) 
and B^ — > K^ip, suggest that one needs ( 8 - 25 % ) non-factorization con- 
tribution to make the model formfactors consistent with data. However, the 
third graph, namely that for the process B~^ K~^ip{2S), suggests a rel- 
atively higher non-factorization but with poorer precision. Concerning the 
decay B^ K^ip{2S) only an upper limit on the branching ratio is available, 
hence we have not plotted the corresponding graph. 



3.2 K*^iIj Decay 

Among B Vtp, the process B^ —>■ K*^tp is the one with the best experi- 



mental data. Values for the branching ratio, (1.58 ± 0.28) x 10~^ [|TT|, and 
polarization, (0.74 ±0.07) are available with reasonable precision for this 
exclusive decay. Hence, we concentrate our attention on this neutral mode. 

As can be seen from (|T^ and ([TB|), non-factorization contributions to this 
decay are parametrized by xai , XA2 ? and xv which correspond to the three 
formfactors Ai,A2, and V. Physics here is more involved. The reason is that 
we have more unknowns than we have constraints. 
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The simplest approach is to assume that only one formfactor has non- 
factorization contribution. We consider this case first. In each of the five 
models, we found that allowing only one of XAi,XA2)and xv to be nonzero, 
one could fit the branching ratio B{B ^ K*^iIj) with appropriate amounts of 
non-factorized contribution. For example, in BSW I model {xai = 0.071 ± 
0.024, XA, = 0,Xv = 0) or {xa, = 0, Xa, = 0.46 ± 0.05, xv = 0) or {xa, = 
0, xa2 = 0, = 0.46 ± 0.1) would fit the branching ratio. However, if we 
fit the polarization measurement, we find that there are only two possibilities, 
iXA, = 0.086±0.098, XA, = 0, Xv = 0) and {xa, = 0, Xa, = 0.32±0.10, xv = 
0) in BSW I model. There are no solutions to polarization data when {xAi = 
0, XA2 = 0, Xy 7^ 0) as shown in Fig. ^ where we show that there are no 
acceptable values of xv that produce large enough polarization in any of the 
five models we have considered. 

Next, we assume that non-factorization is present in both Ai and A2 form- 
factors while Xv is zero. Then the branching ratio demands that the allowed 
region in xai — XA2 space lie between two ellipses as shown in Fig. ^ for each 
of the five models. The region allowed by the polarization measurement lies, 
as shown in Fig. ^ between two pairs of open curves. Thus, in general, there 
are four solutions where the domain allowed by polarization overlaps with the 
domain allowed by the branching ratio as shown in Fig. |^. 

Yet another way to display our results is shown in Fig. ^ where we have 
plotted the region in x, y plane allowed by polarization data. In the factoriza- 
tion approximation {xai = 0, XA2 = 0, Xv = 0) the allowed region is bounded 
by the two curves marked A. If we allow a 5% non-factorization in xai only 
(i. e. xai = 0.05, XA2 = 0, Xy = 0), the allowed region is bounded by the 
two curves marked B. The region bounded by the curves C corresponds to 
(XAi = 0.1, XA2 = 0, xy = 0). The case (xai = 0.1, XAa = -0.03, xv = 0) 
gives the region bounded by the curves D. The values x and y in the five mod- 
els are also shown in Fig. |. The values of xai , XA2 ? a^id xv are chosen only 
to illustrate the effect of non-factorization on longitudinal polarization. 



4 Decays 

Not much is known about the decays of meson. To the best of our knowl- 
edge, only the longitudinal polarization of B^ (pijj (0.56 ± 0.21) has been 
measured [0], albeit with a large uncertainty. However, the results of the 
previous sections regarding the decays B~^ — ^ K^ip, B^ — >■ K*^ip, . . . can be 
used to study the unknown B^ decays if we assume that the amount of non- 
factorized contribution is approximately independent of the light flavor. 

We start with 5° — > {ri,ri'){ip,ip{2S)) decays. For these, the branching 
ratios as a function of xfi , are plotted in Fig. ^ for each of the five theoretical 
models. In order to get some feel for numbers, the branching ratios averaged 
over the predicted values in the five models are shown in Table |^ for %, 10 % 
and 20 % non-factorization contribution. 
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Regarding the process — * (pip, we show in Fig. |^ the regions in x, y 
space allowed by the longitudinal polarization data. The region bounded by 
contours A corresponds to the factorization approximation while that bounded 
by contours C correspond to 5 % non-factorization contribution in xai- The 
region bounded by curves B corresponds to (xai = 0, XA2 = —0.03, Xv = 0)- 
This figure shows that even though the polarization data has low precision, 
factorization approximation does not work for all of the theoretical models 
considered. 

As we did for the other decays in Table ^, we show in Table § the 
branching ratios averaged over the five models and polarization of the processes 
-Bg (pip and B^ (pip{2S) for different values of non-factorization param- 
eters. In Table ^ we have tabulated the ratios B{Bs (pip)/ B{B^ — ^ K^ip) 
and B{Bs (pip)/B{B^ K*^ip) for different amount of non-factorization 
approximation. In this last table the values of B{Bs (pip) were taken form 
Table ^ wheras B{B^ — >■ K^ip) and B{B^ — > K*^ip) were given the experi- 
mental values PT| . 



5 Discussion 

Due to the fact that only two measurements, B{B K*ip) and Pl{B K*ip) 
are available, one can at best derive a constraint between any two of the three 
no n- factorization parameters xai , XA2 ; and xv- The values of these parameters 
depend on the theoretical values of the formfactors Ai, and V . Due to 
the fact that of the three parameters, a, h and c ( eq.(^) - ( P3D ), a is the 
largest, it is most economical (largest effect for the least amount) to put non- 
factorization in xax- We have shown that (5 - 10)% non-factorization in xax 
allows the model formfactors to be consistent with the branching ratio and 
longitudinal polarization data. 

As for the decay B Kip, the five theoretical models we have chosen 
require anything from 8% to 25% non-factorization contribution depending on 
the model. 

Assuming that the amount of non-factorized contribution in {B^,B^) de- 
cays is independent of the light flavor, we have studied the decays Bg 
{rj, ri')ip, (jj, rj')ip{2S), (pip and (pip{2S). To the best of our knowledge, only the 



longitudinal polarization in B^ (pip has been measured |]13[, though with 
large errors. This data is consistent with factorization in BSW I and BSW 
II models though not with the other three models we have chosen. A higher 
statistics measurement would better test factorization in theoretical models. 

The branching ratios for Bg {ri,ri')ip, Bg {ri,ri')ip{2S), Bg (pip and 
Bg (pip{2S) are also quite sensitive to non-factorization. A measurement of 
these branching ratios will be very welcome in testing our ideas. 

A.N.K thanks the Natural Science and Engineering Research Council of 
Canada for the award of a research grant which partially supported this work. 
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Table 1: Model predictions of formfactor Fi{q^) at = or m'^(2S)- 
CDDFGN model, rj stands for rjs, the octet member. This scheme cannot 
handle rji, the flavor singlet. 





BSW I 


BSW II 


CDDFGN 


AW 


ISGW 


B+ - 




0.565 


0.837 


0.726 


0.542 


0.548 


B+ - 


► K+'il]{2S) 


0.707 


1.31 


0.909 


0.678 


0.760 






0.49 


0.726 


0.771 


0.534 


0.293 






0.613 


1.41 


0.964 


0.668 


0.475 






0.411 


0.609 




1.06 


0.463 


B's^ 


W(2^) 


0.514 


0.954 




1.33 


0.752 



Table 2: Model predictions of Ai{m^) , A2{rnV) , and V{mV) formfactors 





A2 V 


X y 




BSWI 
BSW II 

CDDFGN 

AW 

ISGW 


0.458 0.462 0.548 
0.458 0.645 0.812 
0.279 0.279 0.904 
0.425 0.766 1.19 
0.316 0.631 0.807 


1.01 1.19 
1.41 1.77 

1.00 3.24 
1.80 2.81 
2.00 2.56 


Bl ^ H 


BSWI 
BSW II 

CDDFGN 

AW 

ISGW 


0.374 0.375 0.466 
0.374 0.523 0.691 

0.265 0.279 0.919 
0.449 0.703 1.34 
0.237 0.396 0.558 


1.00 1.24 
1.40 1.85 

1.05 3.47 
1.56 2.98 
1.67 2.35 
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Table 3: Average branching ratios predicted by the theoretical models for three 
choices of xfi- 





Factorization 
xlO-^ 


Xf, = 0.1 

xlO-4 


Xf, = 0.2 

xlO-4 


B{B^l - 




0.31 


1.41 


3.1 


B{B^s - 




0.15 


0.64 


1.5 


B{B's - 


rj'ip) 


0.71 


3.1 


7.1 


B{B's - 


rj'i;{2S)) 


0.27 


1.2 


2.7 



Table 4: Average branching ratios and polarization predicted by the theoretical 
models for the processes — * (/>■?/' and B^ (f)ilj{2S). 





Factorization 


XA, = 0.05 
Xa^ = 


XA, = 0.1 

XA2 = 


XA, = 0.1 
XA, = -0.05 


B's ^ 0^ 


BR 
Pol. 


0.21 X 10-=* 
0.32 


0.54 X 10-3 
0.55 


1.1 X 10-3 
0.65 


1.4 X 10-3 
0.72 


5° ^ 0V'(2^) 


BR 
Pol. 


0.14 X lO"-'* 
0.32 


0.32 X 10-3 
0.46 


0.6 X 10-3 
0.52 


0.67 X 10-3 
0.57 



Table 5: Average branching ratios predicted by the theoretical models for 
the process B^ — > (pijj divided by the central values of the experimental mea- 
surements of B{B'^ K'^tli) and B[B^ K*^'il)) for some chosen values of 
non- factorization parameters. 





Factorization 


XA, = 0.05 


XA, = 0.1 


XA, = 0.1 






XA2 = 


XA2 = 


XA2 = -0.05 




0.21 


0.54 


1.1 


1.4 


B(B+^K+i>) 


BiB'i^.pi,) 


0.13 


0.34 


0.70 


0.89 


B(B°-^K*°tp) 
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Figure 1: Quark flow diagram for tlie two body decay of B meson. 
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Figure 2: Allowed region (bounded by the two curves) of function of 

Fi(m^) defined by — >■ K^ip. The dots show the model predictions of the 
formfactors; left: AW, ISGW, BSW I, CDDFGN, BSW II. 
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Figure 3: Branching ratios as a function of xfi predicted by each modeL 
Horizontal hues define the Branching Ratio to one standard deviation. 




Figure 4: Polarization for the process — >■ K*^ip with {xAi = XA2 = 0), 
plotted as a function of xv for each model. Horizontal lines define the measured 
value to one standard deviation. See Fig. ^ for legend. 
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Figure 5: Regions in xai —XA2 space bounded by experimental data on branch- 
ing ratio (ellipses) and polarization (open pairs of curves) for the process 
— > K*^ip taking {xv = 0). Values of xai and xa2 allowed by both data he 
in the four areas where the region between the ellipses overlaps with the space 
between the two pairs of open curves. 
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Figure 6: Regions in x, y plane allowed by experimental data on polarization 
for the process — > K*^il). Contours A bound the region for the case (xai = 
XA2 =Xv = 0), B:(xai = 0.05, XA2 = Xv = 0), C:(xai = 0.1, XA2 = Xv = 0), 
= 0.1, XA2 = ~"0.03,xy = 0).The dots represent predictions of the 
theoretical models. 
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Figure 7: Branching ratios as a function of xfi predicted by each model. In 
CDDFGN model, r] stands for 773 and there is no prediction for rj' . See Fig. ^ 
for the legend. 



15 







ODDFGk^ 












BS^ 


ISGW \ 
II \ 








\ 


BSW I 










ABC 




A B C 




0.5 


1 


7.5 2 2.5 


3 



X 



Figure 8: Regions in x,y plane allowed by experimental data on polarization 
for the process -Bg (pip. Contours A bound the region for the case {xai = 
XA2 = Xv = 0), B:(xai = 0,xa2 = -0.03, xy = 0), C:{xa^ = 0.05, xa2 = 
Xv — 0).The dots represent predictions of the theoretical models. 
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